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Levy-based growth models 
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In the present paper, we give a condensed review, for the nonspecialist reader, of a new modelling 
framework for spatio-temporal processes, based on Levy theory. We show the potential of the 
approach in stochastic geometry and spatial statistics by studying Levy-based growth modelling 
of planar objects. The growth models considered are spatio-temporal stochastic processes on the 
circle. As a by product, flexible new models for space-time covariance functions on the circle 
are provided. An application of the Levy-based growth models to tumour growth is discussed. 
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1. Introduction 

Stochastic spatio-temporal modelling is of great importance in a variety of disciplines 
of natural science, including biology (Cantalapiedra et al. (2001), Brix and Chadocuf 
(2002), Fewster (2003), Gratzer et al. (2004)), image analysis (Feideropoulou and Pesquet-Pop 
(2004)), geophysics (Calder (1986), Lovejoy et al. (1992), Sornette and Ouillon (2005)) 
and turbulence (Schmiegel et al. (2004), Schmicgcl et al. (2005)), to name just a few. In 
particular, the modelling of tumour growth dynamics has been a very active research area 
in recent years (Delsanto et al. (2000), Peirolo and Scalerandi (2004), Pang and Tzeng 
(2004), Schmiegel (2006)). In most of the above-cited works, the model is given implicitly 
and the resulting dynamics are difficult to control explicitly. However, for applications 
and for the theoretical understanding of the modelling framework being employed it is 
essential to connect the ingredients of the model with dynamical and spatial proper- 
ties of the system under consideration. Furthermore, for a parsimonious description of 
systems which arc different with respect to the dynamics and physical mechanisms un- 
derlying the dynamics, it is desirable to have access to a flexible and, at the same time, 
mathematically tractable modelling framework. 

Levy-based models provide a promising modelling framework to meet these require- 
ments concerning flexibility and dynamical control. Until now, Levy-based models have 
mainly been used for describing turbulent flows (Barndorff-Nielsen and Schmiegel (2004), 
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Schmiegel et al. (2004), Schmiegel et al. (2005)). In the present paper, we show that Levy- 
based spatio-temporal modelling also has important applications in stochastic geometry 
and spatial statistics. The main focus is on Levy-based growth models, but we will also 
briefly touch upon another emerging area of application, that is, Levy-driven Cox pro- 
cesses. We expect that the Levy-based approach will have many applications in stochastic 
geometry and spatial statistics. 

Levy-based spatio-temporal models arc constructed from Levy bases, that is, infinitely 
divisible and independently scattered random measures. This terminology was recently 
introduced in Barndorff-Nielsen and Schmiegel (2004). One example of such a model 
describes the growth of a planar star-shaped object, using its radial function Rt(4>) at 
time t and angle 4>. Here, Rt{4>) is the distance from a reference point z to the boundary 
of the object at time t in the direction <j> 6 [— 7t, 7t); see also Figure 2 below. The time 
derivative of the radial function is of the form 



where Z is a Levy basis, A t {4>) C [— tc, 7t) x (— oo, t] is a subset of the past of time t, a so- 
called ambit set (cf. Barndorff-Nielsen and Schmiegel (2004)), ft{-',4>) is a deterministic 
weight function and fit a deterministic function. (In Latin, ambitus means either (1) the 
bounds or limits of a place or district, or (2) a sphere of action, expression or influence.) 
The induced model for the radial function is of the same type. An important advantage 
of these models is that explicit expressions for 



can be derived in terms of the components of the model. This part of the paper is a natural 
continuation of the work initiated in Schmiegel (2006), which was mainly directed toward 
an audience of physicists. An introduction to different growth modelling approaches, 
including a short treatment of the Levy-based approach adopted in the present paper, 
may be found in the SemStat contribution Jensen et al. (2006). 

The paper is organized as follows. Section 2 provides some background on Levy bases, 
which is the essential component of the modelling approach. In Section 3, Levy-based 
spatio-temporal models arc reviewed, while Levy-based growth models are studied in 
Section 4. Section 5 contains explicit results for the covariance functions. In Section 6, 
an application of the Levy-based growth models to tumour growth is discussed, while 
problems for future research are collected in Section 7. 

2. Levy bases 

This section provides a brief overview of the general theory of Levy bases, in particular, 
the theory of integration with respect to a Levy basis. For a more detailed exposition, see 
Barndorff-Nielsen and Schmiegel (2004). As mentioned in the Introduction, a Levy basis 
is an infinitely divisible and independently scattered random measure. Comprehensive 
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accounts of the theory of independently scattered random measures may be found in 
Kallenberg (1989), Rajput and Rosinski (1989) and Kwapien and Woyczynski (1992). 

Let 1Z be a Borel subset of R d , let B = B(1Z) be the Borel subsets of 1Z and let Bb = 
Bb(lZ) denote the class of bounded elements of B. A collection of random variables Z = 
{Z(A) : A G B} or Z = {Z{A) : A G Bb} is said to be an independently scattered random 
measure if, for every sequence {^4 n } of disjoint sets in B (resp. Bb), the random variables 
Z(A n ), n = l,2, are independent and Z({J A n ) = J2 Z(A n ) a.s., where, in the case 
Z = {Z(A) : A G Bb}, we further require 1J A n G Bb- We need to distinguish between the 
two cases B and Bb because the sums ^ Z(A n ) must also be controlled in case where Z 
can take both positive and negative values. If, moreover, Z(A) is infinitely divisible for 
all A G B or Bb, Z is called a Levy basis (cf. Barndorff-Nielsen and Schmiegel (2004)). 

For a random variable X, let us denote the cumulant function logE(e lAX ) by C(A|X). 
When Z is a Levy basis, the cumulant function of Z{A) can, by the famous Levy- 
Khinchinc representation, be written as 

C{\%Z{A))=i\a{A)-} 1 \ 2 b(A)+ /(e iAu -l-iAul hlil] (u))i7(du,A), (1) 

where a is a signed measure on B or Bb, b is a positive measure on B or Bb and U(du, A) 
is a Levy measure on R for fixed A and a measure on B or Bb for fixed du. The measure 
U will be referred to as the generalized Levy measure. The Levy basis Z is said to have 
characteristics (a,b,U). If b = 0, then L is called a Levy jump basis and if U = 0, then 
L is a Gaussian basis; see the examples below. It follows from (1) that any Levy basis 
Z can be expressed as the sum of a Levy jump basis Z\ and an independent zero-mean 
Gaussian basis Z2. 

Without loss of generality (for details, see Rajput and Rosinski (1989)), we can assume 
that there exists a measure fj, such that the generalized Levy measure factorizes as 

tf(du,d£) = V(du,£)M(d£). 

where V(du,^) is a Levy measure for fixed £. Furthermore, the measures a and b are 
absolutely continuous with respect to the measure [i, that is, 

o(d£) = a(0Md£), HdO = 

and a and 6 are uniformly bounded by some constant C > 0. One possible choice (but 
not the only one; see Rajput and Rosinski (1989)) for /j, is 

n(A) = \a\(A) + b(A)+ [ (lAr 2 )U(dr,A), 
Jr 

where \a\ denotes the total variation measure of a and A denotes minimum. 
Let Z'(^) be a random variable with the cumulant function 

C(AtZ'(0)=iA5(£)-iA 2 fe(£)+ f (e iA "-l-iAttl [ _ lil] (u))^(dtt,0- 

Jr 
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factorizable Levy basis and Z'(£) = Z' also does not depend on £. If, moreover, \i is 
proportional to the Lebesgue measure, then Z is called a homogeneous Levy basis and 
all finite-dimensional distributions of Z arc translation invariant. 

The usefulness of the above definitions becomes clear in connection with the integration 
of a measurable function / on 1Z with respect to a Levy basis Z . For simplicity, we denote 
this integral by / • Z . Important for many calculations is the following equation for the 
cumulant function of the stochastic integral / • Z (subject to minor regularity conditions, 
cf., for instance, Barndorff-Nielsen and Thorbj0rnsen (2003)) 



The result (3) can be heuristically derived from (2). A similar result holds for the loga- 
rithm of the Laplace transform of / • Z (assumed to be finite), 



The function K will, in the following, be called the kumulant function. 
We will now give a few examples of Levy bases. 

Example 1 (Gaussian Levy basis). If Z is a Levy basis with Z(A) ~ N(a(A),b(A)) 1 
where a is a signed measure and & is a positive measure, we call Z a Gaussian Levy basis. 
The Gaussian Levy basis has characteristics (a, b, 0) and its cumulant function is 



C(X t Z{A)) = i\a(A) - \\ 2 b{A). 
We have Z'{£) ~ JV(o(0, &(£))> that is, C(\ % Z'(0) = iA «(0 - |A 2 &(£). Furthermore, 



Note that / • Z ~ N(f ■ a,f 2 ■ b). 

Example 2 (Levy jump basis). A Levy basis is called a Levy jump basis if the char- 
acteristics of the basis are (a,0,U). In Table 1, we specify the functions V and a for 
three important examples of Levy jump bases: the Poisson basis, the Gamma basis and 
the inverse Gaussian basis. We also list the distribution of the random variable Z'{t»), its 
cumulant function, mean and variance. All parameters are positive. 

Note that if Z is a Poisson basis, then Z(A) ~ Po(/i(A)) with probability function 







C{\\f-Z)= j C(A/(0 \ Z'(0)/i(dO - iA(/ • a) - \\ 2 (f ■ b) 



(5) 



e 



^li(A) 



X 
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Table 1. The definition of three Levy jump bases - the Poisson basis, the Gamma basis and 
the inverse Gaussian basis - and the distribution of Z'(Q, with the corresponding cumulant 
function, mean and variance 

Poisson Gamma Inverse Gaussian 



V(du,0 


Si{du) 


a«) 


1 


z'(Z) 


Po(l) 


c(\tz'(0) 


c iA -I 


E(Z'(0) 


1 




1 



1 R+ (u)0u- 1 e- a ® u du ^§=1 R+ (u)u-^ 2 e^ 1 ^-' 2 ^ u dz, 

r(/J,a(0) ig(v,i(0) 

-/31og(l-^j) ,7(0(1-^1^) 



P 

7(0 



If Z is a Gamma basis with a(£) = a, then Z(A) ~ r(/3/z(A), a) with density 



r(^(A)) 



/9M(A)-l e -a» ) x > Q 5 



while if Z is an inverse Gaussian basis with 7(£) = 7, then Z(A) ~ IG(t?/u(A),7) with 
density 



W (^)e w(A)7 _ 3/2 f 1,, ,^2 _i 2 , 
4= x 3/2 cxpl~-(( w (A)) 2 x l +-y 2 x) 



x>0. 



The Poisson, Gamma and inverse Gaussian Levy bases are examples of the random 
G-measures introduced in Brix (1999). These measures are purely discrete and can be 
written as (the Levy-Poo representation) 



Z(A)=a (A)+ / xN(dx,A), (6) 
where N is a Poisson basis on M. + x 1Z with intensity measure U and 

a {A) = a(A)- [ xU(dx,A). 



Note that equation (6) can also be written as 



Z(A) = a (A) + £ (7) 
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where $ is a Poisson point process on R + x 1Z with intensity function U . If / is a 
measurable function on 7Z, we have 

f-Z = f-a + u /(0- 

(u,{)€4 

Finally, it should be noted that any Levy process {Zt}teR induces a Levy basis Z on 
ft by 

Z((a,b])=Z b -Z a , a,beR. 

3. Levy-based spatio-temporal modelling 

Let us consider a random variable X t (a) depending on time t and a position <r in space. 
In the following, we will assume that (<r,t) £ 7Z = S x R, where iS C R n , say. A Levy-based 
spatio-temporal model for X — {Xt(cr) : (cr,t) £ 72.} is based on the intuitive picture of an 
ambit set A t (cr) associated with each point (o~,t) £ 1Z, which defines the dependency on 
the past at time t and position a. The ambit set A t (o~) will always satisfy the following 
conditions: 

{a,t)eA t (a), 
A t (a) C5x (-oo,t]. 

An illustration is shown in Figure 1. The linear spatio-temporal Levy model for X = 
{Xt(a) : (cr,f) £ 1Z} is then defined as 

X t (cr)= [ /t(6<r)Z(d£), (8) 

where Z is a Levy basis and ft(£', o~) is a deterministic weight function, which is assumed 
to be suitable for the integral to exist. The process 

X = {exp(X t (o-)):(<T,t)£ll} 

is said to follow an exponential spatio-temporal Levy model. 

The spatio-temporal Levy models can be viewed as generalizations of the familiar 
moving average processes in time series, extended (i) from discrete to continuous time, 
(ii) from one dimension (time) to space-time, preserving a notion of causality in time 
(i.e., the ambit set lies in the past of time t) and (iii) from Gaussian white noise to more 
general infinitely divisible processes. 

Using the key relation (3) , we can derive expressions for moments in the linear spatio- 
temporal Levy model. Thus, we find 



E(X t (a))= [ f t (&o-)E(Z' (£))(,((%) 

JA t {a) 



(9) 
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and 

V(*t(a))= / / t 2 (^ ;( 7)V(Z'(0)^(dO, (10) 
where V is the notation used for variance. The covariances are of the form 

Cov(X tl (cr 1 ) ) X t2 (<7 2 ))= f M^^ft^^mZ'iOMdO- (11) 

J At 1 (a- 1 )nAt 2 (iT2) 

If the weight function is constant, that is /t(£;c) = /, and if the Levy basis Z is factor- 
izable, then (11) reduces to 

Cov(X tl (a 1 ),X t2 (a 2 )) = fY(Z')^(A tl (a 1 )nA t2 (a 2 )). (12) 

In this case, the covariance structure depends only on the u-measurc of the intersection 
of the two ambit sets. 

Equation (4) enables us to calculate arbitrary mixed nth order moments of X t (a) = 
exp(X t (a)). If the moments are finite, then 

E(A (*i) X tn K)) = exp K ^ f t . (£; a, -)l At . M (Q t Z'($J M (d£)J ■ (13) 

The corresponding expressions for the mixed nth order moments of X t (a) are obtained 
from 



nX tl M X tB (<r„)) - gx^~dK E{ ^ (fTl) ^ (an)) 



, (14) 

=A„=0 
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where 

(ffi) Xt (a n )) = exp U K [f^ \ 3 f t] (£; a^A, . (f) t Z'(ol /"(^) j • 

(15) 

The relative second-order moments of X t (o~) have a particularly attractive form 

expl / g(£;t 1 ,t 2 ,a 1 ,a 2 )n(d£) ), (16) 

\JAt, (<7l)nA t ,((72) / 



E(X tl (o-i)A- ta (o-a)) 



E(X tl (cri))E(X t2 (cr 2 )) ViA tl (<Ti)nA t2 (<r 2 ) 

where 

s(£;*x,*2,0'i,0'2) 

= A'((/ tl (£; (n) + f t2 (£; <r 2 )) } " (6 *i) t ^'(0) - (6 ^2) t Z'(0)- 

In the simple case where the weight functions arc constant, that is, /t(£;c) = / for all 
(<r,i) <G 1Z and £ € 71, and where the underlying Levy basis is factorisable, Z'(£) = Z', 
(16) reduces to 

exp(e/i(A tl (CTi)nA t3 (cra))), (17) 

where C = K(2f\ Z') — 2K(f $ Z') . For a factorisable Levy basis Z and a constant weight 
function, one can express 



(18) 



E^VO) E(l t A ;K)) 

in terms of different overlaps of the corresponding ambit sets (Schmicget et al. (2005)) 



4. Levy-based growth models 

In this section, we demonstrate the potential of the Levy setup in stochastic geometry 
and spatial statistics by constructing Levy-based stochastic models for growing objects. 
We focus on planar objects, but generalisations to higher dimensions are straightforward. 
We denote the planar object at time t by Y t C M 2 and will assume that Y t is compact 
and star-shaped with respect to a point z G Yt for all t. The boundary of the star-shaped 
object Y t can be determined by its radial function R t = {Rt(4>) : 4> € [— 7t, 7t)}, where 

Rt{4>) = max{r : z + r(cos <f>, sin </>) € Yt}, 06 [— 7t,7t) 

(cf. Figure 2). 
The growth rate will be described by the equation 

^-R t (<t>)=M4>) + f / t (f;0)Z(df). (19) 
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Here, the deterministic function /i t : [—71, n) — ► R contributes to the overall growth pattern 
while the stochastic integral determines the dependence structure in the growth process. 
The ambit set A t {4>) C [— 7t,7t) x (— oo,t] relates to past events, ft{-',4>) '■ I^ 7 ^ 7 *) xM^M 
is a deterministic weight function (assumed to be suitable for the integral to exist) and 
Z is a Levy basis on [— 7t, 7t) x BL The weight functions and ambit sets must be defined 
cyclically in the angle such that the radial function Rt{4>) becomes cyclic. In the following, 
all angular calculations are regarded as cyclic. 

Note that for nonnegative weight functions and nonnegative Levy bases, Jji? t (0) > 0. 
In other cases, equation (19) still has a growth interpretation if the right-hand side of 
(19) is non- negative with probability one. 

Using (9), the mean growth rate becomes 



dt 



M<t>) 



In the special case where Z is a zero-mean Gaussian Levy basis, m(4>) is indeed the 
mean growth rate at time t in direction <f>. In other cases, /^t(</>) must be chosen such that 
the mean growth rate becomes as desired. There is a large literature on deterministic 
modelling of growth. A classical example is the Gompertz growth rate specified by 



:( ^Rt(<j>) 



~(l-exp(-7*)) 

7 



r;exp(— jt) 



(cf., e.g., Steel (1977)). 

The ambit set A t ((j)) plays an important role in this modelling approach and affects 
the degree of dependence on the past. The extent of the dependence on the past may be 
specified by the minimal time-lag T(t) such that 

A t {4>) C [-7T, n)x[t-T(t),t], <j) G [-7T, 7t) . 

For an illustration, see Figure 3. Note that it follows from the fact that Z is indepen- 
dently scattered that the random growth rates at time t\ and t-z are independent if 




Figure 2. The star-shaped object Yt is determined by its radial function Rt(4>) at time t and 
angle 4>. 
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min(ti,i2) < max(ii — T(ti),t2 — Tfo))- The actual form of the ambit set A t (4>) will 
depend on the specific growth process being modelled. A number of examples are given 
below. The induced correlation structure will be discussed in more detail in Section 5. 
A discrete version of (19) with a Gaussian Levy basis has earlier been discussed in 
Jonsdottir and Jensen (2005). 

For the interpretation of (19) as a growth model, it is helpful to represent the ambit 
set as a stochastic subset of the growing object. This is possible if the stochastic time 
transformation t — > -/?*(</>) is non-decreasing for each <p£ [— 7t,7t). We can then represent 
the ambit set A t (4>) subset of Y t , 

A t (0) = { (R s (0) cos 0, R s (0) sin 0) : (0, s) £ A t (0)}. 

It follows from the fact that A t ((f>) C [— n, n) x (— oo,i] that A t (4>) is actually a subset 
of Y t . Furthermore, since (<f>,t) £ A t (4>), the set A t (4>) touches the boundary of Y t at the 
point (Rt(4>) cos <fi, Rt(4>) sin <f>). It is the 'events' in A t ((f>) that influence the growth rate 
at time t in direction <j>. Figure 4 illustrates the set A t (</>). 

In the particular case where Z is a Poisson basis and ^ the associated Poisson point 
process on [— 7t, 7t) x K, we can represent that part of the spatio-temporal point process 

, arrived before time t, 

*t = {(9 h t i ):t i <t}, 

as a subset of Y t : 

*t = {(R ti (0i) cos °h R U (Si) sm9i) : U < t}. 

We can think of ^>t as consisting of locations of outbursts at time points before t. Finally, 
if we let 

ft ( (s cos 9, s sin 0) ; 4>) = f t ( (6», s) ; </>) , 




4> tf>' 



Figure 3. Two ambit sets A t (<f>) and A t /(<j)'). Note the cyclic definition in the angle. The 
vertical lines are <j> = — n and <f> = rc, respectively. 
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then the fundamental equation (19) can be written as 

J^tfo) £ / t (l;0). (20) 

|e* t nA t (0) 

According to (20) , the growth rate at time t in the direction </> depends on the outbursts 
at time points before t which lie in the stochastic neighbourhood At (</>). This Poisson 
model is closely related to other recently suggested growth models (cf. Section 7.1). 
Under (19), the induced model for Rt{<f) will be of the same linear form since 

R t {4>) = M<i>) + M<i>)+ f I f s {^)z{di)As 

JO JA 3 (<t>) 

= R a {<p) + nM)+ f U^4>)zm, (21) 

JA t (<f>) 

where Rq is the radial function at time t = 0, 

Jk(<t>) = / fi>s(4>)ds, 
Jo 

A t (4>)= (J AM 

0<s<t 

and 

= [ l AsW (0f s fc<t>)ds. (22) 
Jo 




Figure 4. Illustration of the stochastic representation At(<f>) (shown hatched) of the ambit set 
A t (</>) = {{9, s):\e-<t>\<Q,t-T<s<t}. 
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Note that the ambit sets associated with the radial function itself are increasing, that is, 

ti<t 2 ^A tl {(t))<zA t2 {4>). 
If t — T(t) is a non-decreasing function of t, then, because Z is independently scattered, 
R< tl -T( tl ) = {Rt 2 (4>) :t 2 <h- T(ti), 4> 6 [-71, tt)} 

and 

R >tl -R tl ={R t2 (<t>)- RtA4>)-h >t u <j>& [-tt,7t)}, 

are independent. 

The representation (21) is, of course, not unique. If, in particular, 

A t (ct>)=B t nC4 > , (23) 

then 

A t {4>) = B t nC 4> , 

where 

0<s<t 

and we may choose, instead of (22), 

/t(£0)= / i B5 (e)/ s (C;0)ds. 
Jo 

Note that the only difference between the two versions of /t(£; 4>) is the indicator function 
lev 

In some cases, it might be more natural to formulate the model in terms of the time 
derivative of ln(i? t ((/>)), 

!(m(i^))) = ^)+ / A(6$£(d£). 

In this case, the induced model is an exponential spatio-temporal Levy model, 

R t {4>) = Ro{4>)e^{M<f)+ I M^4>)z(do). 

\ JA t {4>) / 

Mixed moments of Rt{4>) can be derived using the results in Section 3. 

The choices of Levy basis Z, ambit sets A t (4>), weight functions /*(£;</>) and /J,t(4>) 
completely determine the growth dynamics. These four ingredients can be chosen arbi- 
trarily and independently, which results in a great variety of different growth dynamics. 
We will now give a number of examples. 



74 K.Y. Jonsdottir, J. Schmiegel and E.B. Vedel Jensen 

Example 3. Consider a Levy growth model for the time derivative of the radial function 

^ t R t {4>) = Z{A t {<t>)\ (24) 

where Z is a Poisson Levy basis with intensity measure concentrated on [—71, 7t) x R + of 
the form 

g(s)dsM, £=(9,s). 
Note that the corresponding point process in the Euclidean plane 

{(scos9,ss'm9):(9,s) is a support point of Z} 

constitutes a Poisson point process with intensity measure 

IMI 

In particular, if g(s) = as, a > 0, then the Poisson point process in the plane is homoge- 
neous. 

Let the ambit sets be given by 

A t (<j>) = 1(9, s) : \9 - (j>\ < ® , max(0, t-T)<s<t 

Represented as subsets of the Euclidean plane, they will, as t — oo, approach rectangles 
of side lengths 28 and T. Note that we can write the ambit set as 

A t {^) = B t nC (t> , 

where 

B t = { (9, s) : max(0, t - T) < s < t} 

and 

C <j> = U9,s):s>-,\9-(l>\<-\u{{9,s):0<s<- 
The mean growth rate at time t and in the direction cj) is, for t > T + ||, 

li{A t {4>)) = 2Q I 9 -^ds. 

Jt-T s 

If g(s) = as, a > 0, then the mean growth rate is constant. Figure 5 shows simulations of 
this model with constant mean growth rate. 
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Example 4- The size of the ambit sets plays an important role in the control of the 
local and global fluctuations of the boundary of the object Y t . As an example, let us 
consider a Levy growth model of the form 

Rt(<l>)=fH + Z(A t {<f>)), (25) 

where 

A t (4>) = {(0, s) : \6 - ft < &{s),t - T(t) <s<t}. 
In Figure 6, simulations are shown under this model, using a normal Levy basis with 

Z(A)~N(0,o- 2 ll(A)) 

and n equal to the Lebcsguc measure on 1Z. Note that (j,(At(<f>)) does not depend on 4>. 
The simulations are based on a discretization of Z on a grid with At = 1 and A</> = . 
The upper and lower row of Figure 6 show simulations for two choices of angular extension 
of the ambit set at three different time points. The angular extension of the ambit set is 
9(s) = for the upper row and 9(s) = ^ for the lower row. For the smaller angular 
extension, we observe localized fluctuations of the profiles, but the global appearance is 
circular. For the larger angular extension, the fluctuations are on a much larger scale and 
the global appearance is more variable. 



Example 5. In this example, we study a model similar to the one described in Exam- 
ple 4, but now with a Gamma Levy basis. The model equation is 

R t {4,)=JX t + Z{A t {4>)), (26) 

where At((f>) is defined as in Example 4, 

Z(A)~T(J3ii(A),a), 




Figure 5. Simulation of the Levy growth model (24) for the derivative of the radial function at 
time points t = 75, 100 and 125, using a Poisson Levy basis. The parameters of the simulation 
are g(s) = 10s, T = 1 and 6 = 1/2. 
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Figure 6. Simulation of the Levy growth model (25) at time points t = 20, 45 and 80, using a 
Gaussian Levy basis. The upper row and lower row show simulations of two choices of the angular 
extension of the ambit set O(s) = and G(s) = respectively. Otherwise, the parameters of 
the simulation are H20 = 16, /Ms = 24, fiso = 32, a 2 — 1 and T(t) = t/5. 



and fi denotes Lcbcsguc measure on 1Z. The parameters a, (3 and fit are chosen such 
that E(i?t (</))) and Y(Rt (</>)) are the same as in the previous example. Here, we have 
used Table 1, together with (9) and (10). Accordingly, the parameters are chosen such 
that 



jit = fit - (Ty/f3fJ,(A t (0)), 



The only free parameter is j3 > 0, which determines the skewness of the Gamma distri- 
bution of Z(At((f>)). For large values of j3, the distribution will resemble the Gaussian 
distribution. 

The resulting simulations for /3 = 1 are shown in the upper and lower rows of Figure 7 
for two choices of angular extension of the ambit set, 6(s) = ^ and 9(s) = j, respec- 
tively. Note that more sudden outbursts are seen compared to the previous example. 



Example 6. In Figure 8, we show simulations from the Levy growth model 

R t (4>) = f(4>)(»t + z(A t (m, 



(27) 
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I = 20 




I. = 45 




/ = 80 




I ■ : 20 




I = 45 




I = 80 




Figure 7. Simulation of the Levy growth model (26) at time points t = 20, 45 and 80, using 
a Gamma Levy basis. The upper row and lower row show simulations of two choices of the 
angular extension of the ambit set 0(s) = j^i an d ®( s ) — § ; respectively. Otherwise, (3 = 1 and 
the remaining parameters are determined by the parameters used in Example 4. 



- 20 




I - 15 




.' - 80 




Figure 8. Simulation of the model (27) at time points t — 20, 45 and 80, using a Gaus- 
sian Levy basis, with parameters as specified in Example 4. The weight function is given by 
/ t (0) = O.35exp(^). 



where fit , A t (</>) and Z are specified as in Example 4 and 

f t {4>) = 0.35 exp 



— 7t 



71 



Clearly, the growth of the object is asymmetric. The weight function /t(</>) puts more 
weight on the angle (f>o =0. 
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5. The induced covariance structure 



Levy-based growth models lead to flexible new models for space-time covariance functions 
on the circle, as we shall see in this section. The results presented here are of general 
interest for spatio-temporal processes on the circle. 

We will derive expressions for Cov(i? fl (</>i), Rt 2 (<fc)) under various assumptions on the 
Levy basis Z, the ambit sets A t ((f>) and the weight functions /{(£; (j>). We will concentrate 
on the Levy growth model (21) of linear form for R t . Since we now are interested in 
covariances, it suffices to look at the model equation 

Rt{<t>)=l mA)z{&0, 

JA t {4>) 

where we, for simplicity, have omitted the bar on the ambit set and weight function. The 
covariance structure of Rt(4>) is then given by (cf. (11)) 

Cov(R tl (<j )1 ) 7 R t2 (fo))= [ / tl (C;0i)/ t2 (C;02)V(z'(e)HdO- (28) 

JA tl (4> 1 )nAt 2 {4>2) 
Throughout this section, we will assume that 

A t (<i>) = (ct>,0)+A t (0), 
/ t (£;0) = /t((|0-0|,s);O), (29) 

for all £ = (0, s) £ TZ and (<j>,t) £1Z. These conditions ensure that Cov(Rt 1 ((f>i),Rt 2 (<p2)) 
only depends on the cyclic difference between 4>i and 2 . Accordingly, the spatio-temporal 
process 

{iJt(0):t€R,0G[-7t,7t)} 

will be second-order stationary in the space coordinate, but not necessarily in the time 
coordinate. 

We will first consider the case where the angular extension of the ambit set is the full 
angular space, but the weight functions are arbitrary. Second, we consider the case of 
constant weight functions, but quite arbitrary ambit sets. 

5.1. Ambit sets with full angular range 

In this subsection, we consider ambit sets of the form 

A t (<f ) ) = [-n,n)x [t-T(t),t]. 
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In order to express the formulae as compactly as possible, we use, in the proposition 
below, the notation t\ fl t 2 for the time points shared by A tl {-) and A t2 {-), that is, 



\ otherwise, 



f 

where 

ii =max(*i — T{t\),t 2 — T(t 2 )) and t 2 = min(ii, t 2 ). 

Using this notation, we can derive the following convenient and general expression for 
the covariances. 

Proposition 7. Let us assume that the ambit set is of the form A t ((j>) = [— 7t, n) x [t — 
T(t),t] for all (<f>,t) G K and let 

OO 

/ t (f; <f>) = al(s) + J2 a k(s) C ™(H0 ~ 0)), (30) 
fe=i 

£ = (0,s), be the Fourier expansion of the weight function. The spatio-temporal covari- 
ances are then 

00 

Cov(R tl (0i), R t2 (0 2 )) = 2T (t 1 ,t 2 ) + J2 Tfc(*i. *2) coB(ft(& - 0a)), (31) 

fe=i 

T k (t l ,t 2 )=n [ a%(s)a%(s)g(s)ds. 

Proof. The proof is straightforward. First, note that the actual form (30) of the Fourier 
expansion of the weight function is a consequence of (29). We obtain that 

Cov(R tl (^),R t2 (<f> 2 )) 
•/A tI (^i)nA t2 (&,) 

= 7t 2 f < ( S )a* 2 ( S ) 5 (s) ds + J2( [ ( S )4 2 (s) 5 (s) d.s) cos(fc(0 1 - 2 )) . 

Jtint 2 k=l \Jtint 2 J 

□ 

Exam/pie 8. Suppose that the weight function is of the form (30) with a.j.(s) = if k ^ 1. 
Then, 

Cov(i? tl (0 1 ),i? t2 (0 2 )) = 7Tcos(0 1 -0 2 ) f a t 1 1 {s)a\ 2 (s)g(s)ds, 

Jtxnt-2 
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Since the covariance is a product of a spatial term and a temporal term, this model is 
separable (cf. Stein (2005) and references therein). The sign of the covariance may be 
positive or negative. 

Note that, according to (31), the covariance Cov(i? tl (<j>i), Rt 2 (<fo)) depends on <j>i and 
cf>2 only via \<j>i — <j>2\ - For some choices of model parameters, the covariance also becomes 
stationary in the time coordinate. For instance, if g(s) = 1, T(t) = T and a k (s) — b k (t — s), 
then we have 

^T+mm(ti-t 2 ,0) 

Tk(ti,t 2 ) = TC b k (u)b k (t 2 -ti+u)&u. 

J max(ii — 12 ,0) 

The induced model (31) for the covariance function is not, in general, separable in 
the sense that the covariance function can be written as a product of a term depending 
only on t\ and t 2 and a term depending only on (jii and </>2- This may be regarded 
as a strength of the model because separable covariance functions are often believed 
to give too simplistic a description of spatio-temporal data (cf., e.g., Stein (2005)). If, 
nevertheless, such simplifying assumptions are made, we obtain the following results. 

Corollary 9. Let the assumptions be as in Proposition 7. Assume that a k (s) =a\- The 
spatial correlations are then determined by the Fourier coefficients of the weight function 

fR U \ b U V - Cov(Wi),W2)) _ 2K) 2 + E^iK) 2 cos(fc(0i"02)) 
P (R t (M,R t (h)) ■ A(Wi))v(W2)) - 2^ + lXrK) 2 ' 

//, in addition, a\ = btc k , then the covariance model (31) is separable. Furthermore, the 
spatial correlations p(R t {<f>i) , Rt{4>2)) do not depend on t, while the temporal correlations 
are determined by T(t) and the function g: 

L ™ 9( s )ds 

p{R tl R t , ((b)) - — - t Unt2 t . 

LC-TCtx) 9(s)ds ■ j t 2 2 _ T{t2) ff(s)ds] 1 /2 

The covariance model (31) provides a possibility for extending stationary covariance 
functions on the circle (spatial covariance functions) to a spatio-temporal context. When 
Rt is a stationary process on the circle, its covariance function can be expressed as 

00 

CovtRtfai), RM) = 2Xl + J2 X k C0S (H<t>i ~ <t>2))- (32) 

fe=i 

Such a covariance function can be obtained by choosing the Fourier coefficients of the 
weight function to be 



400 = 4 = -j= 



1/2 

Al/ / g(s)ds 

t-T(t) 
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Note that there is still freedom in the modelling by choosing an arbitrary time-lag T(t) 
and function g. 

Example 10. The pth-order model for a stationary covariance function on the circle, 
described in Hobolth et al. (2003), has 

Aq = A* = 0, A| = [a t + l3 t {k 2 P - 2 2 p)}-\ k = 2,3, . . . . 

The model is called a pth-order model because it can be derived as a limit of discrete pth- 
order Markov models defined on a finite, systematic set of angles (cf. Hobolth and Jensen 
(2000)). This covariance structure is obtained by choosing 



a*(s) = a\(s) =0, 

r t 1-1/2 

4(*)= " 



7t / #(s)ds 

Jt-T(t) 



[a t +p t (k 2 P~2 2 P)}-^ 2 , fc = 2,3,. 



If a t and /3 t are proportional, then the simplifying assumptions of Corollary 9 are fulfilled. 
In Jonsdottir and Jensen (2005), this model has been used for the time derivative of the 
radial function. Only Gaussian Levy bases are considered and neighbouring time points 
are assumed to be so far apart that the increments can be regarded as independent. The 
more general approach of the present paper allows for temporal correlations. Under the 
assumption a\ — btc^, the temporal correlations are particularly simple. For instance, 
suppose that T(t) = 1 and t 2 — 1 < ti < ^2- We then get for g(s) = ae~ bs , a, b > 0, 

p{R tl (0), R t2 = -^[eCVWi-tO+i- _ e (-i/2) b ( h -t 2 ) h 
while, for g(s) = as a , a > 0, a > 1, 

ff +1 - (t 2 - 1) Q+1 



P {RtMM<t>))- 



Only in the first case are the temporal correlations always stationary. 
5.2. Constant weight functions 

In this subsection, we consider the case of constant weight functions. Without loss of 
generality, we assume that ft(£,',4>) = 1 an d so (28) reduces to 



Cov{RtM,R t M) = / V(Z'(0)fi(dO=V(Z')(,(A tl (<t>i)nAtM), 
J A tl (<t>i)nA t2 (<t> 2 ) 

where the last equality holds if the Levy basis is factorisable. 
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It is not difficult (but sometimes tedious) to find explicit expressions for Cov(i? tl (</>i), 
Rt 2 {4>2)) for specific choices of ambit sets. One simplifying assumption is to focus on 
ambit sets of the form 



Usually, it is easier to find expressions for the temporal covariances than for the spatial 
covariances. 

Evidently, (33) implies that Cov(i? tl (</>i), Rt 2 (^2)) > 0, which may be a severe restric- 
tion for the spatial covariances. In the proposition below, the spatial covariances are 
expressed in terms of the function delimiting the ambit set. The proposition provides 
insight into the class of spatial covariances that can be modelled using this approach. 

Proposition 11. Let /i(d£) = g(s)dsd9 for £ = (9,s). Let us suppose that there exists a 
continuous function h t : [—71, 7t) — > M with the properties 



A t ( ( j ) ) = B t nC <P7 



where 



Bt 



{(9, s) : max(0, t ~ T(t)) <s<t}, 
{(M):|^-0|<e(*)}. 



ht(6) = h t (-9), 

h t is decreasing on [0,7t), 

h t (0)=t 



(34) 



such that 



A t (0) = {{e,s):ht{n)<s<ht(0)} 



(cf. Figure 9). Let 




If the Fourier expansion of h t is (ht(4>) = ht{—4>) ) 



00 




(35) 



fe=0 



then 




(36) 



where 





Proof. It is not difficult to show that 

/ — 7T+(/>/2 />7T 
ht(0)d9 + 2 ht(0)d0-2nht(n), 0e[O,7t). (37) 
-71 ^0/2 

Using (35), we find that 

- 4 E y sin (4) +2TC E ^- 27t E ^> if«^e[0,7r] 

k odd fe odd even 

4 ^ i sin ^ +27t ^ 7 *-27t £ 7 *, if 0£ [-71,0]. 

A: odd A: odd k even 

The result is now obtained by deriving a Fourier expansion of the latter expression and 
comparing with (36). □ 

Example 12. In the particular case where g(s) — 1 and 

h t (</>) = h t ((f)) = 7q + 7* cos <f>, 
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we find that 



2n- 



16 

71 

16 1 



7i 



-7i, 



i = 1,2,. 



It follows that 



(A5) _1 =at + Ai s 



1.2, 



(38) 



where at = —Tij (I67') and /3t = n/ (47'). Under the assumption of a normal Levy basis, 
(38) is a special case of the pth order model considered in Jonsdottir and Jensen (2005) 
with p = 1 and a proportional to [3. Note that requirements (34) imply that 7q = t — 
and 7* > 0. It does not seem to be possible to obtain pth order models with p > 1 using 
this approach. 



6. An application to tumour growth 

In Schmiegel (2006), snapshots of a growing brain tumour in vitro were analyzed using 
the approach described in this paper; see Figure 10. The data were first studied in Bru 
et al. (1998). 

A detailed initial analysis of the covariance structure showed negative spatial covari- 
anccs and a need for modelling both small and large scale fluctuations in the growth 
process. The model used was an exponential spatio-temporal Levy model of the form 

{/•t-io(i) pit 
(H + a(t) / / cos(0-<9)Z(d.sd<9) 

Jt-T(t) J-n 

(39) 

pt p<p+h t {s-t+t Q (t)) ~| 

+ (3(t) / / Z(dsd6) \. 

Jt-t a {t) J 4>-ht(s-t+t (t)) J 




Figure 10. Profiles of a growing brain tumour in vitro at time points t = 21, 25 and 55. 
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Table 2. The estimated parameters for model (39), using a Gaussian Levy basis 
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t 


T(t) 


to{t) 


«(*) 


m 




21 


21 


1!) 


0.04 


-0.033 


0.19 


25 


25 


17 


0.02 


-0.033 


0.19 


55 


18 


4 


0.01 


-0.067 


0.23 



Here, ht is a deterministic and monotonically decreasing function defined on [0,io(*)]j 
satisfying h t (to(t)) = and /i*(0) = <j>o(t)/2. Accordingly, the weight function is of the 
form 

ft(€;<t>) =a(i)cos(^-6»)l [t _ T(t)! t_ to( t )] (s)+/3(<)l [ t_ to(t)jt] (s)l [0 j lt(s _ t+ t o(t))] (|(/)-e'|). 

The associated ambit set is shown in Figure 11. In Schmiegel (2006), a Gaussian Levy 
basis was used and the function h t was assumed to be of the form 



h t (s) 



2 2t (ty 



3£[0Mt)]- 



The parameters of model (39) were estimated by the method of moments, using the results 
given in Section 3. The estimated parameters arc given in Table 2 and a simulation under 
the model with a Gaussian Levy basis is shown in Figure 12. 

Here, we will study the use of Gamma and inverse Gaussian Levy bases. Simulations 
under the latter basis are shown in Figure 13. The inverse Gaussian Levy basis is chosen 
such that E(R t (4>)) and Y(Rt(4>)) are the same as in the case where a Gaussian basis is 
used. The upper row of Figure 13 shows simulations where r\ = 316 and the lower row 



i 



t - Ht) 



i-T(i) 




4>-n<t>-<l>o{t)/2 d> <t> + M* )/' 2 <t> + * 
Figure 11. The ambit set At((j>) for the model defined by (39). 
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shows a simulation where 77 = 5. For n = 316, the inverse Gaussian Levy basis provides 
fits of a similar quality as the normal basis, but for r) = 5, more outbursts are observed, 
as is the case for the data. The difference is due to the fact that the inverse Gaussian 
distribution has heavier right tails for the latter choice of parameters. It should be noted 
that all of the profiles simulated under model (39) using the Levy basis mentioned in 
this section show slightly more fluctuations on a local scale than the observed profiles. At 
present, we do not know whether this feature is caused by non-perfect model selection 
and estimation of parameters or artefacts due to the discretization in the simulation 
procedure. 

7. Discussion 

In the present paper, we have given a condensed review of Levy-based spatio-temporal 
modelling and shown its potential use in stochastic geometry and spatial statistics by 
developing Levy-based growth models and space-time covariances on the circle. Below, 
we discuss further perspectives and topics for future research. 

7.1. Related growth models 

In the growth literature, there is a variety of growth models for objects in discrete space 
(cf., e.g., Bramson and Griffeath (1981), Qi et al. (1993), Lee and Cowan (1994), Kansal 
et al. (2000) and references therein). An important early example is the Richardson 
model, introduced in Richardson (1973). Here, the growth is described by a Markov 
process. For a growing object in the plane, the state at time t is a random subset Y t of 
Z 2 consisting of the 'infected sites'. An uninfected site is transferred to an infected site 
with a rate proportional to the number of infected nearest neighbours. It can be shown 
that if Yo consists of a single site, then Y t /t has a non-random shape as t — ► 00. Note that 
the growth model described in (20) of the present paper may be regarded as a continuous 
analogue of the Richardson model. 



t = 21 

o 



Figure 12. Simulation of model (39) for time points t = 21, 25 and 55, using a Gaussian Levy 
basis. 
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I = 21 

o 



f =21 

o 



Figure 13. Simulations of model (39) for time points t — 21, 25 and 55 using an inverse Gaussian 
Levy basis with 77 = 316 (upper row) and 77 = 5 (lower row). 



A related growth model in continuous space has recently been discussed in Deijfcn 
(2003). For planar objects, the model is constructed from a spatio-temporal Poisson 
point process on R 3 , 

* = {(ii,ti)}. 

The random growing object Y t C M. 2 is a subset of 

IJ B( Xi ,r), 

{i: ti<t} 

constructed such that Y t is always connected. Here, B(x,r) is a circular disc with centre 
x and radius r. In this model, U is thought of as a time point of outburst and Xi as 
the location of the outburst in the tumour, say. A closely related discrete-time Markov 
growth model has been discussed in detail in Cressie and Hulting (1992). This model can 
be characterized as a sequence of Boolean models, 

Y t+1 =U{B(x i ,r):x i eY t }, 

where {xi} is a homogeneous Poisson point process in M 2 ; see also Cressie and Laslett 
(1987) and Cressie (1991a, 1991b). 

An issue of interest in growth modelling is the asymptotic shape of the growing object 
(cf., e.g., Durrett and Liggett (1981) and Deijfen (2003)). It is expected that it is also 
possible to obtain asymptotic results for Levy-based growth models using the fact that 
Levy bases are independently scattered random measures. 
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7.2. Levy-driven Cox processes 

Another interesting application of spatio-temporal Levy models in spatial statistics is 
provided by Levy-driven Cox processes (cf. Hellmund (2005), Hellmund et al. (2007) and 
Prokesova et al. (2006)). As an example, we may use exp(X t (cr)), where X t (a) is given 
in (8), as driving random intensity for a spatio-temporal Cox process. If Z is a Gaussian 
Levy basis, then the resulting Cox process is log-Gaussian. Another example is obtained 
by assuming that Z is a positive Levy basis, the Cox process driven by X t (a) then being 
a spatio-temporal shot-noise Cox process. Cox processes have been studied intensively in 
recent years (cf. Brix (1998), Wolpert and Ickstadt (1998), Brix (1999), Brix and Diggle 
(2001), Brix and M0ller (2001), Brix and Chadoeuf (2002), M0ller (2003)). 



7.3. Estimation of model parameters 

It still remains to develop inference procedures for Levy-based growth models (and, more 
generally, for Levy-based spatio-temporal models). There are several interesting problems 
here, including nonparamctric estimation of the ambit sets. 

A Fourier expansion of the radial function may be useful when making inference about 
the shape of the growing object (cf., e.g., Alt (1999) and Jonsdottir and Jensen (2005)). 
Let us consider the Fourier coefficients of Rt(4>), 

1 f n 1 f n 

Ai = - / R t U) cos(fc0) d(j), Bl = - / R t {4>) sin(fc0) d<f>, 

k = 0, 1, . . . . Under the assumptions of Proposition 7, it can be shown that 

\ al(s)cos(k6)Z(d6ds), B\= / / a' fc (s)sin(fc0)Z(d0da), 

-nJt-T(t) J-7iJt-T(t) 

so the Fourier coefficients also follow a linear spatio-temporal Levy model. It can be 
shown that, for k ^ j, t, t' > 0, 

Cov(Al , A*' ) = Cov(B* , Sj' ) = Cov(4 , Bj ) = 

and 

Cov(4 , At ) = Cov(i?* , Bt ) = r k (t, t% 

where Tk(t,t') is given in Proposition 7. 

In the case where Z is a Gaussian Levy basis, this means that {j4i}tgR and {-B|.}igR, 
fc = 0, 1,..., are independent Gaussian stochastic processes with covariance functions 
Tk(t, t'). If one observes A\ and B\ for some time points t = ti, . . . ,t n and some orders 
fc = 1, . . . ,K t , the likelihood function is tractable. 
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